An exact quantification of backreaction in relativistic cosmology 
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An important open question in cosmology is the degree to which the Friedmann-Lemaitre- 
Robertson- Walker (FLRW) solutions of Einstein's equations are able to model the large-scale be- 
haviour of the locally inhomogeneous observable universe. We investigate this problem by consid- 
ering a range of exact n-body solutions of Einstein's constraint equations. These solutions contain 
discrete masses, and so allow arbitrarily large density contrasts to be modelled. We restrict our 
study to regularly arranged distributions of masses in topological 3-spheres. This has the benefit of 
allowing straightforward comparisons to be made with FLRW solutions, as both spacetimes admit 
a discrete group of symmetries. It also provides a time-symmetric hypersurface at the moment of 
maximum expansion that allows the constraint equations to be solved exactly. We find that when 
all the mass in the universe is condensed into a small number of objects (< 10) then the amount of 
backreaction in dust models can be large, with 0(1) deviations from the predictions of the corre- 
sponding FLRW solutions. When the number of masses is large (> 100), however, then our measures 
of backreaction become small (< 1%). This result does not rely on any averaging procedures, which 
are notoriously hard to define uniquely in general relativity, and so provides (to the best of our 
knowledge) the first exact and unambiguous demonstration of backreaction in general relativistic 
cosmological modelling. Discrete models such as these can therefore be used as laboratories to test 
ideas about backreaction that could be applied in more complicated and realistic settings. 

PACS numbers: 98.80.Jk,04.20.Jb 



I. INTRODUCTION 



The visible matter in the universe contains structure on 
a wide range of scales, from stars and planets (~ 10^°m) 
to voids and super-clusters (^ lO^^m), (see e.g. [T]). In 
standard general relativistic cosmological modelling the 
problem of dealing vifith this complexity is often circum- 
vented by assuming that the stress-energy tensor of an 
inhomogeneous universe can be approximated by an 'av- 
eraged' stress-energy tensor, usually assumed to be rep- 
resentable by a simple perfect fluid, while leaving the left 
hand side of Einstein's equations unchanged. The so- 
lution to these equations is then taken to represent the 
geometry of spacetime on the scale over which the av- 
eraging was performed. If this scale is large enough to 
result in a homogeneous and isotropic spacetime, then 
we will refer to this geometry as the "background" . 

The great benefit of this approach is the enormous sim- 
plification it entails, by allowing cosmology to be done on 
any scale, without having to worry about the enormously 
complicated task of taking into account the wide range 
of structures that exist in the universe. The drawback is 
that it is not always clear to what extent the expansion 
of space in the background geometry is representative of 
the large-scale expansion of space in the actual universe. 
For a given model, we will refer to the difference between 
these two things as "backreaction" , which can be consid- 
ered as the influence of inhomogeneities at different scales 
on large-scale cosmological dynamics. As large-scale ex- 
pansion can be related to observables such as luminosity 
distance and redshift [2H1|, it is therefore the case that 
an understanding of backreaction could prove essential 



in achieving a fundamental concordance between theory 
and observations [3 [6] . 

This task becomes increasingly important given that 
observational cosmology is now starting to be able to 
make high-precision measurements of a variety of dif- 
ferent astrophysical probes, including the Cosmic Mi- 
crowave Background (CMB) [315], Type la supernovae 
[91110], and Baryon Acoustic Oscillations (BAO) [TT ] IT2 ] . 
When interpreted within the background of the standard 
cosmological model, these very different probes all sug- 
gest that the universe is at present undergoing a phase 
of late-time accelerated expansion. In a genuinely homo- 
geneous and isotropic general relativistic universe filled 
with a perfect fiuid, such behaviour is only possible if 
p < — p/3, with data currently being compatible with 
the presence of a cosmological constant with p = —p. 
The presence of such a fiuid is problematic for a variety 
of reasons, requiring a fine-tuning of perhaps as much as 
1 part in 10^^*^. This immediately raises the question of 
whether the model we are using to interpret the data is 
adequate, a question which needs to be quantified in any 
case if we are to start performing "precision cosmology" 

m 

One particular aspect of this complicated problem 
is the degree to which dust (a fiuid without any self- 
interactions) is appropriate for describing the matter con- 
tent of the universe [14] . This approximation treats grav- 
itationally bound, extended objects (such as galaxies, 
and clusters of galaxies) as a medium that is parametrised 
by only a handful of continuous, and usually slowly vary- 
ing, quantities. Within this framework, the spacetime 
through which photons travel has a fundamentally differ- 



ent type of curvature to that of the real universe [T5] . It 
also ignores the gravitational consequences of the binding 
energies between and within astrophysical systems [IB] . 
A dust description should therefore be considered as an 
approximation only, and must itself be investigated in 
order to ascertain the extent to which it is justifiable. 

We study the problems outlined above by consider- 
ing cosmological solutions that are explicitly composed 
of regularly arranged discrete masses. These solutions 
consist of vacuum everywhere exterior to the sources un- 
der consideration, and are constructed on a manifold 
that has spacelike hypersurfaces that are topological 3- 
spheres. As a result, they admit a group of discrete 
symmetries corresponding to rotations of the 3-sphere 
that leave the positions of the masses invariant. These 
models are among the simplest possible inhomogeneous 
spacetimes that could be considered to exhibit large-scale 
homogeneity and isotropy; They introduce only one ex- 
tra scale into the problem (the inter-particle separation) , 
they are vacuum solutions of the field equations (up to 
the singularities), and they contain within them a time- 
symmetric spatial hypersurface (which greatly simplifies 
the constraint equations). This makes them ideal candi- 
dates for studying backreaction. 

Cosmological solutions of this type have been consid- 
ered before, starting with the seminal work of Lindquist 
and Wheeler [T7]. The approach these authors took was 
to construct a gravitational analogue of the highly suc- 
cessful Wigner-Seitz construction. The basic idea was 
to build lattices from regular cells, and attach the cell 
boundaries tangentially to a background 3-sphere. Dy- 
namics were then imposed by assuming that observers 
that are equidistant between neighbouring sources should 
be falling freely with respect to those sources. This led 
to a dynamical model of the universe in which the masses 
are regularly arranged, and that is taken to be an approx- 
imate solution to Einstein's equations. The Lindquist- 
Wheeler approach has recently been extended to include 
models with spatially flat and open topologies [T5], as 
well as to models with a cosmological constant [19]. The 
optical properties of spatially flat models of this type 
have also been considered in p!5H20| . as have those of 
negatively curved models in pi] . Other attempts to con- 
struct cosmological models for similar conflgurations of 
discrete masses have been made using perturbative ex- 
pansions [U, and Regge calculus [23]. 

While providing insights into the problem, all of the 
models mentioned above involve the use of approximate 
solutions. Our aim here is to take an alternative approach 
to addressing this problem that is exact. This approach is 
based on work initiated by Misner about flve decades ago 
[24j . and that has subsequently been used extensively in 
the study of black hole physics i25-rS7] . The underlying 
idea here is again drawn from analogy with the study of 
electromagnetism, where intuitions about more general 
situations are gained by flrst studying the simple case of 
static configurations of charge. In gravitational settings, 
however, it is very difficult to arrange for configurations 



of isolated masses to be static for a finite interval of time^ . 
In contrast, configurations that are instantaneously static 
occur frequently. The study of such configurations, that 
occur in many systems of physical interest, is referred to 
as geomctrostatics^ . 

In the case of cosmology, which is our main interest 
here, an important example is provided by the instanta- 
neously static (and therefore time-symmetric) hypersur- 
faces that occur at the moment of maximum expansion 
in dust-dominated FLRW solutions with positive spatial 
curvature. Our goal is to construct a set of instanta- 
neously static configurations that are exact, and that 
contain discrete matter sources, rather than continuous 
fluids. We do this by constructing regular tessellations of 
a 3-sphere by using different sets of identical polyhedra, 
each with an identical mass at their centre. We show that 
the method of gcomctrostatics can be extended to solve 
the constraint equations exactly in these cases. Our aim 
is then to use these solutions to study backreaction in 
dust filled cosmological models. We do this by compar- 
ing the time-symmetric spatial hypersurfaces of the exact 
discrete solutions to the corresponding FLRW solutions 
with the same proper mass^, and considering, in partic- 
ular, the extent to which the FLRW solutions emerge as 
the number of sources is increased. This provides a pre- 
cise indication of the degree to which FLRW solutions 
can be used to represent initial configurations for space- 
times that contain discrete objects. A fuller investigation 
of backreaction along these lines would also include con- 
sidering the dynamical evolution of models with discrete 
sources. We leave this for a future publication. 

The plan of the paper is as follows. In Section jTT] we 
recap on perfect fluid FLRW cosmology. In Section III 
we provide the constraint equations on time-symmetric 
hypersurfaces for discrete masses on a 3-sphere. We dis- 
cuss the Schwarzschild solution, and then outline how to 
flnd the solution for arbitrarily many discrete masses in 
a closed space. In Section [TV] we construct all of the pos- 
sible solutions that consist of discrete sources regularly 
arranged on a closed lattice. This allows us to explain 
why there are no 2-inass solutions with spherical topology 
to Einstein's equations that admit a time-symmetric hy- 
persurface. In Section [V] we provide evidence that shows 
that the sources in our discrete models are always sep- 
arated by distances that are larger than their horizon 
sizes at the maximum of expansion. In Section |VI| we 
investigate backreaction in dust-dominated cosmological 
models by comparing the scale of our lattices to dust- 



^ Counter-examples include spacetimes in which the cosmological 
constant is non-zero, and the periodic solution found by Korotkin 
and Nicolai .5B'_. 

^ Examples include dynamical systems that possess a recoUapsing 
phase (as happens when a particle is thrown vertically upwards 
from the surface of a massive body, with less than the escape 
velocity). 



As defined in |29l and in Section IIIB 



filled FLRW solutions of Einstein's equations with the 
same total proper mass. We find that 0(1) deviations 
from the results of Friedmann cosmology can occur when 
the number of masses is small (< 10), but that the scale 
of the solutions converge to the Friedmann values when 
there are very many masses (> 100). In Section VII 



we then consider what happens when we include inter- 
particle interaction energies in our definition of mass. We 
find that the consequences of changing the number of 
masses in the lattice, while keeping the total energy of 
the system the same, can considerably change the scale 
of the hypersurface of maximum expansion. This shows 
that although backreaction in dust models we consider 
is generically quite small, the consequences of ignoring 
interaction energies when using the dust approximation 
can hide potentially interesting effects. In Section VIII 
we present our conclusions. 



II. PERFECT FLUID COSMOLOGY 

The spatially homogeneous and isotropic perfect fluid 
FLRW solutions have a geometry given by 



ds' 



^f{t)de + a\t) (dx^ + h\x)dn^) , (1) 



where f{t) is a free function, dfl'^ — dd^ + sin^(6')d0^, 
and h(x) — sinx, X or sinhx for solutions with spatial 
curvature k = +1, or —1, respectively. The scale factor 
a{t) satisfies the constraint equation 



1^ 
/a2 



Stt 



(2) 



where the dot represents a derivative with respect to t, 
and p = p(t) is the energy density of the continuous per- 
fect fluid, which obeys the conservation equation 



p + 3-{p + p)=Q, 
a 



(3) 



where p — p{t) is the pressure. This completely specifies 
the solution, up to constants of integration. 

For the present study we will be primarily interested 
in models with positive spatial curvature (fc — +1), and 
with a pressureless dust source (p = 0). In this case, if 
we take / = a(i), then the solution can be written as 



, , Stt 1 , ,9 



(4) 



where to and po are constants, and where p{t) — po/a^{t). 
The maximum of expansion can then be seen to occur at 
t = to, and the geometry of the hypersurface t = to can 
be seen to be given by 



df^---—{dx^+sm\dn' 

&TTp{to) 



(5) 



which is rigidly specified once p{to), the energy density 
at maximum of expansion, is known. 



We note that the only dust-filled FLRW solutions 
that admit a time-symmetric hypersurface are spatially 
closed, with a spherical topology. Interestingly, the 
boundary conditions for discrete models with a momen- 
tarily static distribution of sources also seem to be incom- 
patible with open topologies (the reason for this is that a 
1/r source cannot live alone on T^, and a regular lattice 
of identical sources on E^ is also not possible [28j ) . This 
shows a certain qualitative concordance between discrete 
and fluid solutions in the interplay between boundary 
conditions and topology. 



III. CONSTRAINT EQUATIONS FOR 
DISCRETE MODELS 



Our aim here is to obtain exact vacuum solutions of 
Einstein's equations corresponding to regular lattices of 
sources that are instantaneously at rest on a topological 
3-sphere. 

In such a setting the relevant equations to solve are 
the Gauss-Codazzi equations: 



n + K^ 



K,jK'^ = 







{K/~S^K).^ = 0, 



(6) 

(7) 



where TZ is the Ricci curvature of the 3-space, Kij is the 
extrinsic curvature of the 3-space in the 4-dimensional 
spacetime, and K = K^^. The indices i, j refer to co- 
ordinates in the 3-space, and the vertical line denotes 
covariant derivative in that space. 

It is well known that if we choose a time coordinate 
that is specified by the normal derivative to this initial 
hypersurface (such that gtf_i = ~<5„*), then the extrinsic 



curvature can be written as Ki. 



Inst ant a- 



neously static hypersurfaces therefore have Kij = 0, and 
the Gauss-Codazzi equations reduce simply to 



n = o. 



(8) 



A key point here is that for any 3-dimensional geometry 
that satisfies the initial constraint ([8]) there is a unique 
4-dimcnsional spacetime that satisfies the full Einstein 
equations. This is the method of gcometrostatics, pre- 
sented by Misner in [M|, and studied for the case of 
asymptotically flat space in p5H?7] . 

Here we are interested in solving Eq. (Is]) in closed 
spaces. In this case we can make the following ansatz 
for the metric of the spatial 3-section 



d?" 



i/j h.ijdx^dx\ 



(9) 



where ^ — 'ip{x^), and where hij is the metric of a 3- 
sphere with constant curvature TZ. The Gauss-Codazzi 
equations are then satisfied if ip obeys the Helmholtz 
equation: 



where V^ is the Laplacian corresponding to hi 



(10) 



A. A Single Schwrarzschild Mass 



B. Many Schwarzschild Masses 



The Schwarzschild solution has been well studied using 
asymptotically flat solutions in geometrostatics, includ- 
ing in the original work of Misner [M] . It is important to 
note that the Schwarzschild solution is still a solution in 
the present case, as the 3-sphere is conformally flat. In 
this case, however, there will be a different radial vari- 
able and a different functional form for the conformal 
factor Tp. Given the linearity of the Helmholtz equation, 



(10 1, it can be seen that multi-source solutions can be 



constructed by linear superposition of the solutions for 
single sources. 

In asymptotically flat, isotropic coordinates, a static 
spacelike slice of the Schwarzschild solution can be writ- 
ten as 



dl' 



1 



m 
2f 



[df^+r^dn'^) 



(11) 



where ttt, is a constant. Performing the coordinate trans- 
formation f = K tail ^, where X is a constant, it can be 



seen that Eq. (Ill becomes 



dV 



4 



2i^sin| 



{dx^ + si-a^xd^'^)- (12) 



This metric is clearly of the same form as Eq. (|9| , and 
solves Eq. (10), with TZ — 6 (for a unit 3-sphere) and 



V' 



1 m 

3S # 2K sin $ 



(13) 



Actually, both functions A{x) ■— 
(cos |)~i satisfy Eq. ^ with 7^ 



(sin I) 
= 6. 



and B{x) 



This means that although the term in Eq. ( 12 ) that 
corresponds directly to the 2m/r source term from 
Eq. (Ill can be seen to be the one that is proportional 



to A(xj, the term that is proportional to B{x) can also 
be treated as a source, although it is located at the an- 
tipode. This becomes obvious if we place the origin of 
the spherical coordinates at the antipode by the trans- 
formation X — > TT — X, as in that case the roles of the two 
terms in Eq. ^V2\ are interchanged (i.e. A{x) -^ B{x) 
and B{x) -^ a(x))- Thus it appears that placing a 
Schwarzschild source on the 3-sphere induces a mirror 
source at the antipode. Furthermore, these two source 
terms are joined at their horizons, as is clear if we set the 
gauge parameter to K = m/2 (in which case the horizon 
at f = m/2 appears at x = 7r/2 in the hyperspherical 
coordinates). 

One can then consider that as x — ^ tt one approaches 
either the asymptotic region where f — > oo, or that one 
approaches another Schwarzschild mass. These two situ- 
ations are geometrically identical, and there is therefore 
no difference between the one mass solutions and the two 
mass solution (as long as the point at x = '"' can be added 
to the manifold). 



From Eq. ( 12 1 we can now infer that there exist multi- 



source solutions on the 3-sphere that take the form 

df = ^'^ {dx^ + sin^ xd9?), (14) 



where ip = '0(Xj ^j 4>) is given by 



N 



Hx,o,4>) = Y. 



-[2U{x.e,, 



(15) 



where the mass parameters fhi are a set of constants and 
fi{x,(^i4>) = sin(xi/2). Here the Xi refer to new coor- 
dinates {xi,0i,4>i) which are obtained by rotating the 
coordinates (x, ^, 0) in Eq. (14) so that the z'th source 
position appears at Xi = 0- 

It follows from the form of Eq. 



( 14 ) that one cannot 



adjust the fhi and the size of the 3-sphere independently. 
For example, scaling the size of the 3-sphere by a constant 
a^ automatically results in scaling the rhi by a factor a. 
This demonstrates that there exists a certain rigidity in 
the configuration, once the fhi are specified. Specifying 
the value of the fhi is therefore crucial for determining 
any measure of distance in our hypersurface. Once they 
are specified, however, then all measures of distance can 
be calculated uniquely. 

When considering the original single source 
Schwarzschild solution, the parameters fhi {fhi and 
7712 in that case) are equal to the standard Schwarzschild 
mass, fhi = 7712 = ?ti (using the gauge K ~ m/2). 
However, when multiple sources are present, the inter- 
pretation of the fhi parameters will instead be that of 
an effective mass which includes the binding energies 
with respect to all the other objects in the universe 
(cf. [li]). However, their actual, locally measured, 
proper mass'^ must be determined by the requirement 



that the geometry reduces to that given by Eq. ( 12 ) in 
the limit x ~> 0, in analogy with the asymptotically flat 
situation analysed in Ref. [5S]. It can be shown that 
the proper mass defined in this way is equivalent to the 
mass in the Schwarzschild region at the other end of the 
Einstein-Rosen bridge of the black hole. 

To summarise, these two mass definitions are: 

1 The Effective Mass, {mi)cs ^ is defined to be equal 
to the mass parameter of the corresponding source, 
(mi)eff := fill- 



* This term is used in analogy with the proper mass defined by 
Wald |29l p. 126]. In the context of asymptotically flat space- 
times, Brill and Lindquist 25 use the term bare mass, as it is 
used in particle physics. We prefer not to use this term here 
since the situations in cosmology and particle physics are op- 
posites with respect to how masses are measured. In particle 
physics, the observer is outside the system and cannot measure 
the bare mass directly. In cosmology, on the other hand, the ob- 
server is inside the system and therefore measures the bare mass 
(in particle physics terminology). 



2 The Proper Mass, rrii, is defined in such a way that 
the geometry of the space in the vicinity of each 
object is given by Eq. (12 1, in the limit x ^ 0. 



Strictly speaking, one requires a region to exist that is 
infinitely far away from all of the masses in order to have 
an operational definition of the gravitational effect of the 
inter-particle potentials that are included in the effective 
mass. This is clearly impossible in the case of a closed 
space, but this does not stop us from drawing an anal- 
ogy with cases in which asymptotically distant regions 
do exist. 



IV. LATTICE SOLUTIONS 

Using the above ingredients we can now construct a set 
of models with different numbers of regularly arranged 
discrete masses on a 3-sphere. To do this we proceed 
in analogy with the approach taken by Lindquist and 
Wheeler [T7] , so that our models are constructed by con- 
sidering all of the possible regular tessellations of the 3- 
sphere. We then place identical masses at the centre of 
each cell, with the result that the distribution of masses is 
such that the distance between any mass and its nearest 
neighbours is the same for each of them. 

There are six possible regular tessellations of the 3- 
sphere, those with 5, 8, 16, 24, 120 and 600 polyhedra 
|30j . In addition, one can also cover a hypersphere with 
2 balls [3T]. These possibilities are displayed in Table 
|Ij Using these tessellations, we will now proceed to con- 
struct a sequence of exact discrete regular lattice models 
of the universe, with increasing numbers of sources, on a 
time-symmetric 3-sphere. Our discrete models will con- 
sist of n equal masses (where n = 2,5, 8, 16, 24, 120, 600). 
By positioning these masses at the center of each of the 
cells we will then achieve a regular distribution on the 
3-sphere. 



Lattice 


Cell 


Number of 


Structure 


Shape 


Cells 


- 


Ball 


2 


{333} 


Tetrahedron 


5 


{433} 


Cube 


8 


{334} 


Tetrahedron 


16 


{343} 


Octahedron 


24 


{533} 


Dodecahedron 


120 


{335} 


Tetrahedron 


600 



A. The 2-Cell model 

The 2-cell model, consisting of two balls with their 
boundaries identified, is the exceptional case of the struc- 
tures listed in Table |I] Its cells are not regular polyhedra 
(and hence it has no Schlafli symbol), and it was not con- 
sidered in either [17] or [30]. Nevertheless, a lattice with 
2 cells seems like a perfectly legitimate object to con- 
sider^ and has indeed been studied recently in [3T]. The 
authors of this paper found that the geometry inside each 
cell could only be matched at the junction between cells 
if that junction was a horizon. They also found that no 
solutions exist unless the cosmological constant A 7^ 0. 

Let us now consider this 2-cell model using the initial 
value formalism discussed above. The source functions in 
this case are given by /i = sin(x/2) and /2 = cos(x/2), 
with rhi = ifi2 ~ rn. Substituting into Eq. ( |l4| then 
gives 



V- 



2 cos ^ 



2 sin ^ 



(16) 



which is clearly just the Schwarzschild solution given in 
■m/2. As explained above, this model 



121, withX 



Eq, 

can be considered to consist of two sources centred at 
opposite poles and matched (analytically continued in 
fact) across their horizons. Within this interpretation 
the exterior parts of the sources are missing, and there 
are consequently no static regions outside of the black 
holes. Because of the absence of an exterior region be- 
tween the sources, the 2-cell model is unsuitable for cos- 
mology. This result provides an alternate illustration of 
the findings in Ref. [5T] . 



B. The 5-Cell model 

Let us now consider the 5-cell model (or hyperpyra- 
mid), that consists of 5 tetrahedra. This structure, which 
is in fact a 4-simplex, is the tessellation in Table |T] with 
the fewest number of cells that are regular polyhedra. To 
obtain the coordinates of the masses at the centres of the 
tetrahedra one can consider a hyperpyramid in the em- 
bedding space E^, and place a unit 3-sphere inside the 
hyperpyramid such that the two structures touch at the 
centre of each of the 5 tetrahedral cells. Alternatively, 
one could consider the positions of the vertices of the 
dual lattice^, which in this case is another hyperpyra- 
mid. The coordinates of the 5 masses that result are 
given in the Table |ll| below. 



Table I: All possible regular tessellations of the 3-sphere. 
The 'Lattice Structure' is given by the Schlafli symbols {pqr}, 
where p is the number of edges to a face, q is the number of 
faces that meet at a vertex, and r are the number of cells that 
meet at an edge {30f . 



^ This object does, in fact, have a considerable advantage over 
the lattices in Table llj as the geometry around each mass must 
be spherically symmetric, and so is considerably easier to solve 
for. This means that the approximations used by Lindquist and 
Wheeler in 1171 are not required in this case at all. 

® Dual lattices have the vertices and centres of each cell transposed 
with each other. 



To define the model we use spherical polar coordinates 
which are related to Cartesian coordinates of £"* by 



w = cos X 
X = sin X cos 9 
y — sin X sin 9 cos ( 
z = sin X sin 9 sin (j 



(17) 



The source functions, fi, from Eq. (14 1, are then 
f, =sin[icos~i(/ij)] , 



(18) 



where the functions, hi, defined in Eq. (181 are given by 



hi 
h2 = 

h. = 



hi — 



cosx 



15 „ . COSY 

— cos 9 sm Y 



sm X sm 9 sm < 

5 . . cosx 

sm X cos f 



sm X sm 9 sm 



sm X cos t 



cosx 



■ sm X sm 9 sm 



(^-?) 



5 . cosx 

sm X cos f 

The geometry of this model can be visualised by consid- 
ering a slice through it. To this end, consider the surface 
X = Xo for xo — cos^^(— 1/4) w 1.82. Its metric is given 

by 



16 



]dn^ 



(19) 



To get a rough idea of the shape of this hypersurface 
we can think of 9 and (j) as polar angles in E^ and plot 



the surface tp{xo,9,(f)) = p where p^ = x'^ 



y 



This surface is displayed in Fig. la and goes through 4 
of the 5 masses. It should be noted that the surface does 
not represent an isometric embedding of the geometry in 
Eq. (19). However, in regions where the derivatives of ^ 



are small, it does give an approximate representation of 
that geometry. This approximation is therefore best in 
regions which are far from the sources. 



C. The 8-Cell model 

We now proceed in a similar manner to find the geom- 
etry of the 8-cell (or tesseract) model. In this case the 
primitive cell of our lattice is a cube. To find the position 
of the centre of each cell we can again embed the struc- 
ture in £"*, together with a unit 3-sphere. Alternatively, 



Point 


{w, X, y, z) 


(X, 9, <t>) 


(ii) 


(1,0,0,0) 

(-i4^.o,o) 




(o.f.f) 

(cos-^(-|),0,f) 


{Hi) 


(-i-\/i.\/i.o) 




(cos-i(-i),cos-i(-|),0) 


(ill) 


/ 1 Air AT 

l^ 4 ' V "*** ' V 24 ' \ 


/i) 


(cos-i(-i),cos-i(-|),f) 


(v) 


I^^4'^y48'^y24'^ 


Vi) 


(cos-i(-l),cos-i(-|),f) 



Table II: Coordinates (w, x, y, z) of the 5 masses in the embed- 
ding space iJ*, as well as (x, d, 4>) on the background 3-sphere. 
In this table, and throughout, cos~^ refers to the inverse co- 
sine, and not its reciprocal. 



the positions of the masses can be found using the dual 
lattice, which in this case is the 16-cell. The positions of 
the masses are then given as in Table |TTTj 
The fi from Eq. ( 14 ) are then found to be 



/i = sin 

/2 = cos 

/a = sin 

fi = cos 

/5 = sin 

/e = cos 

fr = sin 

/s = cos 



cos ^ (cos sin x) 



- cos ^ (cos 9 sin x) 



cos^ ^ (cos (f> sin 9 sin x) 
cos" ^ (cos (j) sin 9 sin x) 
cos" ^ (sin sin 9 sin x) 
cos" (sin (j) sin 9 sin x) 



The first two of these functions are identical to those 
corresponding to the 2-cell model, considered above. 

We can again visualise the geometry of this model by 
considering a slice though it. In Fig. [lb] we show the 
surface with x — ^/2, which passes through 6 of the 8 
masses. The radial position of the surface at any given 
(6', (p) is proportional to Tp. 



D. Models with 16-600 Equally Spaced Masses 

We can construct the other discrete models, made us- 
ing 16, 24, 120 and 600 equally spaced masses, by pro- 
ceeding in a similar way to the cases discussed in detail 
above. We shall not present the details of these construc- 
tions here, but to help with their visualisation we display 
slices through these structures in Figures Tc][Tf It can be 
seen that as the number of masses increases, the shape 
of each of these structures becomes increasingly spher- 
ical, while the tubes become thinner. This corresponds 
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metric on the 2-space. This condition is satisfied if 



Table III: Coordinates {w,x,y,z) of the 8 masses in the em- 
bedding space _£*, as well as {x,d,(t>) m the lattice. 



4sin2xsine'(l + xf^){l + xfe) {x\4>'^\4> + Xiei'ie - V',x) 
+ip sine [I + X\4> + X\e) X\4>4' + (^ + xf^j X\ee 
+V'cosxsin0[(l + X|0)X|e + ^'^(l + X|0 + X|e) sin0 

-sinx(2 + 3xf^ + (3 + 4xf^)xfe)] 
= -'0sinxcos6'x|e(l + xfe), (20) 



to the spacetime approaching homogeneity as the number 
of masses is increased. We wih use the resuhs obtained 
from studying these larger lattice models in the sections 
that follow. 



LOCATION OF THE HORIZONS 



where for compactness we have introduced the notation 



X\e 



,Xj 



and 



X\4> 



^X,tj>- The positions of 



sm X ' >^ sm X sm f 

the horizons are then approximated by the solution to 
this equation. 



To qualify as cosmological solutions we require that the 
discrete models considered here avoid having any overlap 
in the horizons corresponding to different masses. We 
therefore need to investigate the positions of the hori- 
zons in the models discussed in the previous section, in 
order to see if this criterion is met. In these models the 
location of the event horizon can be approximated by 
marginally trapped surfaces |26j (such surfaces give the 
exact locations of the event horizons if the spacetime is 
static) . 

To find these trapped surfaces let us consider a surface 
given by some function x — x(^)^)- This surface has 
geometry 



In practise Eq. (20) is not easy to solve, but we can 



obtain approximate solutions by looking for minimal 
surfaces of constant Xi for which Eq. (20 1 reduces to 
('0*sin^x).x ~ 0- We expect this to be a good approx- 
imation for all the models considered here, and for its 
accuracy to increase as the number of masses in the lat- 
tice is increased. To establish this result we calculate the 
area of the horizon of each mass using our approximation, 
Aijiin , as well as the horizon area of a Schwarzschild black 
hole with an equal proper mass, ^s- We then calculate 
the ratio (^min — ^s)/^Si which we have displayed in Ta- 



ble IV As expected, the difference is small, and decreases 
as the number of masses is increased. 



da^ 



^\x% + sm'xW + IxfiX^^i^HBd^ 



+^\x\ 



'xsm^0)dcf>^, 



and the unit normal to it is 



"0^ sin X sin ( 



sin'^X sin^6' -I- sin^6'x „ + X 



(1,-X,e,-X,0): 



where commas in subscripts denote partial differentia- 
tion. The extrinsic curvature of this surface is then 



tJ.v 



''IJ,\vi 



where the covariant derivative here is 
with respect to the metric of the 3-space. Transform- 
ing to coordinates a, 6,... on the 2-space, using Kab = 



dx'' dx 



dx^ 

K - 



-Q^K^i,, the trapped surfaces are those that have 
- j'^^Kab = 0, where 7"'' is the contravariant induced 



To check that the horizons are not overlapping we com- 
pare Xmin with half the separation between neighbouring 
sources. Ax. This is also displayed for each of our 6 dis- 
crete models in Table |IV[ It can be seen that Xmin is 
always less than half Ax, and that it decreases as the 
number of masses is increased. Together with the small 
values of (Amm — As)Ms this provides a good indica- 
tion that the horizons of the masses in our models do not 
intersect at the maximum of expansion, thus ensuring 
that our discrete models satisfy a necessary condition to 
qualify as cosmological models^. 



'^ We leave aside for now the more complicated question of whether 
there are additional horizons that could encompass two or more 
masses. We note only that this would appear to be unlikely, 
given that in every case the masses are separated by multiple 
horizon distances. 





(a) A slice through the 5-ceU solution 



(b) A slice through the 8-cell solution 





(c) A slice through the 16-cell solution 



(d) A slice through the 24-cell solution 



0G 



(e) A slice through the 120-ccll solution 



(f) A slice through the 600-cell solution 



Figure 1: Slices through the hypersurfaces of the discrete lattice solutions. The distance from the centre is proportional to t/j. 
The tubes correspond to the locations of the masses. These objects become more and more spherical as the number of sources 
in the lattice increases. 



Number of 

Masses 


Ah 


Xmin 

Ax/2 


5 

8 

16 

24 

120 

600 


0.00114 
-0.00126 

0.000208 
-0.000272 
-0.000749 

0.000601 


0.428 

0.268 

0.173 

0.110 

0.0330 

0.0147 



its maximum of expansion by 



Table IV: The fractional difference between our estimate of 
the horizon size, Amin, and the horizon size of a Schwarzschild 
black hole with the same proper mass, As- Also displayed is 
our estimate of the fraction of the distance to the point half- 
way between masses that the horizon reaches, for each of our 
discrete models. 



VI. BACKREACTION IN DUST MODELS 

We shall now employ our exact discrete solutions in 
order to study backreaction in dust-filled cosmological 
models. This will be done by comparing the scale of 
discrete and continuous models on their time-symmetric 
hypersurfaces, at the maximum of expansion. 

Here we will make use of the notion of 'proper mass' 
in order to compare our lattice models with the dust- 
filled (k = +1) FLRW solutions. That is, for a given 
lattice we will calculate the total proper mass of all the 
sources in the lattice, and we will then compare this lat- 
tice to an FLRW solution with the same total proper 
mass. The motivation for this procedure is that in both 
cases, the proper mass corresponds to locally measured 
masses. Specifically, for the FLRW models, the proper 
mass for a given region can be defined as the integral of 
the energy density over that region. For the lattice mod- 
els, the mass within each cell is identified with the proper 
mass of the source it contains, so that the total proper 
mass of the lattice is given by the sum of the proper 
masses of all the sources within it. 

The geometry of the time- symmetric hypersurface of 
maximum expansion in (fc = -1-1) FLRW solutions is now 
given by Eq. ([5|. On the other hand, the metric cor- 
responding to the 3-sphere containing n discrete masses 
of equal size is given by Eq. (14), with rhi+i = nn for 
The value of m, 



every i < n. I'he value of rrii can then be related to the 
proper mass of each of the objects using the procedure 



outlined in Section IIIB In order to compare continuous 
and discrete solutions we then only need to make sure 
that the mass of the continuous solutions (defined as the 
constant M = pV where V is the total spatial volume of 
the universe) is the same as the sum of the proper masses 
in the discrete solutions. 

In the discrete solutions, consisting of n objects each 
with proper mass m, the total mass is clearly just M = 
nm. In the FLRW solutions the volume of a spatial sec- 
tion of constant t is given by y = 27T^a^, where a — a{t) 
is the scale factor. Thus the energy density for a (fc = +1) 
FLRW solution with the same total mass (M) is given at 



p{to 



M 



M 



27r2a3(to)' 



(21) 



Recalling that in a (fc = +1) FLRW solution the maxi- 
mum of expansion occurs when a? — 3/(87rp), the line- 
element ([5]) can be written as 



dl^ 



16Af2 
97r2 



{dx^ +sin^Xd(^^ 



sin^Xsin^^d^^) . (22) 



We can now compare this geometry with the correspond- 
ing discrete geometry given by Eq. ( 14 ) . 



To proceed we require a measure of the global scale for 
both the discrete and the continuous solutions. In the 
continuous case it is clear what this measure should be, 
as there is only one scale in the geometry (the curvature 
of the 3-sphere). For the discrete solutions, on the other 
hand, the length of a curve of fixed angle on the 3-sphere 
will depend on its particular position, as the geometry of 
the space is inhomogeneous in these solutions. We must 
therefore proceed with some care. Here, for the discrete 
solutions, we propose two possible measures of the 'size' 
of the space. These are: 

Dl The line-element, dl, of a curve at a vertex of the 
lattice. These are clearly distinguished positions 
within the lattice, corresponding to the points that 
are furthest from all masses. 

D2 The length of the edge of a cell. Again, this is 
clearly a preferred curve within the lattice. 

Both of these proposed definitions of size are aimed 
at trying to identify the scale of the lattice structure 
that the masses occupy, as this is the closest thing to 
a "background" that exists in these solutions. Identify- 
ing the size of the lattice was also what was attempted 
by Lindquist and Wheeler in their approximate solution 
with the same configuration of masses [17] . Within this 
context, Definition Dl will turn out to be the most con- 
servative possible comparison of the scales of the discrete 
and continuous solutions, and Definition D2 will be found 
to be not very different. 

Let us first consider Definition Dl. In this case it can 
be seen that the ratio of the line-elements of curves that 
cover the same angle on the 3-sphere in the continuous 
and the discrete solutions is given by 



dl 



discrete 



dFLRw 



16nm l^ 



(23) 



Due to the fact that rhi can be shown to be in direct 
proportion to m, it can be seen that this expression is, in 
fact, independent of m. It is therefore specified uniquely 
by the structure of the lattice (i.e. by n, and the functions 
fi corresponding to the lattice in question). 

The ratio of scales in Eq. (23) is a function of x, & and 



and so varies depending on which point in the discrete 
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Figure 2: The scale of the discrete solutions as a fraction of 
the scale of the continuous solutions, when the total mass in 
the continuous model is taken to be equal to the total proper 
mass of the discrete model. The lower full curve (blue) shows 
the result of using Definition Dl for the scale of the discrete 
solution, and the upper full curve (red) shows the result for 
Definition D2. The dotted curve, which represents the ap- 
proximate Lindquist- Wheeler solution, is shown for compari- 



solution we wish to consider. Using Definitions Dl and 
D2, we find the results displayed in the third and fourth 
columns of Table W\ respectively. Here aQ'^'^''°'° refers 
to the scale of discrete solution, as defined using Dl or 
D2. These results are also displayed graphically in Fig. [2] 
for each of our lattices. It can be seen that the scales 
corresponding to the discrete solutions is always larger 
than that of the corresponding continuous solution, with 
the only exception being the scale of the solution with 
600 masses, when Definition Dl is used. 



0.20 



0.15 



;o.io 



0.05 



0.00 




Figure 3: The ratio of the effective mass, m^ff, to the proper 
mass, m, for each lattice solution. 



Figure ^ shows that a^''''=™*7a^^^w approaches 1 as 
n becomes large, but that this approach is not exactly 
monotonic. In particular, the ratio for the lattice with 
16 masses is a little lower than might have been the case 
for a smooth curve. This behaviour is likely to be a con- 
sequence of the fact that the shape of the cells used in dif- 
ferent tessellations are different. For the lattice solutions 
made from tetrahedra (i.e. those with n — 5,16 or 600 
masses), for example, the approach to the FLRW limit 
does appear to be monotonic. For the sake of comparison, 
we have also shown in Fig.[2]the ratios obtained using the 
approximate solutions of Lindquist and Wheeler |17j . It 
can be seen that the ratios given by the approximate so- 
lutions compare well with those of the exact solutions. 

The outcome of this comparison demonstrates that, 
for large numbers of regularly arranged masses at their 
maximum of expansion, this measure of backreaction in 
dust-fiUcd models is small. 



Cell 


Number of 




I „FLnw ) 

\ '^O / D2 


Shape 


Masses 


Tetrahedron 


5 


1.321 


1.360 


Cube 


8 


1.258 


1.291 


Tetrahedron 


16 


1.061 


1.097 


Octahedron 


24 


1.083 


1.099 


Dodecahedron 


120 


1.033 


1.034 


Tetrahedron 


600 


0.996 


1.002 



Table V: The scale of the discrete solutions as a fraction of 
the scale of the continuous solutions, when the total mass in 
the continuous model is taken to be equal to the total proper 
mass of the discrete model. The two definitions of scale in the 
discrete solutions have both been calculated. 



VII. EFFECTIVE MASSES AND 
INTERACTION ENERGY 

An advantage of the discrete models we consider is 
that they provide us with a framework within which it is 
possible to discuss the role of inter-particle interactions 
in cosmology. Given that these are the interactions that 
govern the structure and evolution of the real universe, 
we consider it to be an important undertaking to obtain a 
deeper understanding of them. In particular, interactions 
are ignored when approximating the universe as being 
filled with dust, and so these considerations potentially 
allows us further insights into the behaviour of spacetimes 
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Figure 4: The scale of lattice solutions with the same total 
effective mass, normalised so that ao = 1 for the 5-cell, and 
using Dl for the scale of the discrete solutions. 



filled with discrete objects. 

Here we will begin by comparing the notions of proper 
mass and effective mass for each of our discrete solutions, 
as defined in Section llll Bl We remind the reader that our 
definition of effective mass is based on analogy with the 
asymptotically flat case studied by Brill and Lindquist 
[25], as it is problematic to define it operationally in the 
closed spherical settings we are currently studying. Now, 
the values of these masses can be set to any given value 
for each lattice, but the ratio of effective mass to proper 
mass must take a particular constant value for any given 
configuration. We display this ratio for each of the six 



possible tessellations in Table VI and show it graphically 
in Fig. [3] It is clear that the ratio of effective mass to 
proper mass increases as the number of masses in the 
lattice is decreased. This is due to decreasing contri- 
butions from the interaction energies between particles, 
which correspondingly increase the value of the effective 
mass of each source. 



Cell 


Number of 


"^cff 


Shape 


Masses 


m 


Tetrahedron 


5 


0.20 


Cube 


8 


0.11 


Tetrahedron 


16 


0.045 


Octahedron 


24 


0.029 


Dodecahedron 


120 


0.0052 


Tetrahedron 


600 


0.0010 



In Fig. [4] we compare the scale of the hypersurface of 
maximum expansion for each of our discrete solutions, 
when the total effective mass of the lattice is kept con- 
stant. As the effective mass includes the inter-particle in- 
teraction energies, this procedure of comparing the scale 
of solutions with the same total effective mass is equiv- 
alent to enforcing the condition that the lattices being 
compared should contain the same total energy (that is, 
the same total proper mass plus interaction energies). 
Fig. |4] then shows that the scale of the hypersurface of 
maximum expansion grows approximately linearly with 
the number of masses, n. Changing the method of deter- 
mining the scale of the discrete solutions from Dl to D2 
has very little effect. 



The results in Fig. [4] can be understood by considering 
Fig. [3] - as the number of masses decreases the ratio of 
the magnitude of interaction energies to proper mass also 
decreases. The interaction energies, however, are nega- 
tive, so if we consider a thought experiment in which we 
deform the 8-cell into the 5-cell we have to increase the 
total amount of energy available for interactions. The 
only source of energy in the solution is the proper mass 
of the particles, and so a lattice with fewer masses must 
have a smaller total proper mass if the total energy in 
the system is to remain unchanged. Smaller total proper 
mass corresponds to a smaller scale for the hypersurface 
of maximum expansion, as was shown in SectionlVl] This 
result suggests that one can substantially increase the 
scale of a closed space by dividing up the mass in that 
space into smaller packets, although we have strictly only 
shown this is the case on time-symmetric hypersurfaces^. 
It should also be noted that this result relies on the defi- 
nition of effective mass that we have generalized from the 
asymptotically flat setting in order to define the interac- 
tion energies. The validity of this approach is difficult 
to confirm operationally in a compact space, and so the 
interpretation of the dramatic increases in scale seen here 
should be treated with care. 



Such increases in scale are entirely absent in dust mod- 
els, as interaction energies are not included in the energy 
budget in that case. This does not, however, invalidate 



any of the conclusions of Section VI It remains the case 
that back-reaction is small in models where interaction 
energies are ignored. What has changed for the lattices 
in Fig. [4] is that the proper mass in each lattice is in- 
creasing as the number of masses is increased, so the 
dust-dominated FLRW solution that one should compare 
with these models is also changed. 



Table VI: The ratio of the effective mass, m^s, to the proper 
mass, m, for each lattice. 



* For a dynamical cosinological model we no longer necessarily 
have energy conservation, and so the situation could be more 
complicated. 
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VIII. CONCLUSIONS 

We have studied the emergence of spatially closed 
Friedmann solutions from inhoniogeneous solutions that 
contain increasing numbers of regularly arranged discrete 
masses in topological 3-spheres. This has been done using 
exact methods, and therefore allows for exact measures 
of backreaction in dust dominated cosmological models. 

More specifically, by considering the instantaneously 
static hypersurfaces at the moment of maximum expan- 
sion we have found that universes that contain only a 
small number of mass concentrations (< 10) can be 10% 
or more larger than the corresponding dust dominated 
FLRW solutions of Einstein's equations. However, for 
universes that contain very many masses (> 100) there 
is very little quantitative difference (< 1%) between the 
scale of maximum of expansion in the discrete and con- 
tinuous solutions. These results are based on compar- 
ing dust dominated continuous models (which ignore 
the inter-particle interaction energies) with momentar- 
ily static discrete models that contain the same total 
'proper mass' (which also ignores the interaction energies 
between the different masses). While backreaction in this 
case is quite small, it may well be that the inter-particle 
interactions will come into play in a more substantive way 
when we allow for dynamics. That is, backreaction could 
still have a big effect when we consider the evolution of 
the universe. 

The results we find are in good keeping with the ap- 
proximate solutions of Lindquist and Wheeler that use a 
gravitational analogue of the Wigner-Seitz construction 
[T7] , as is demonstrated graphically in Fig. [2] They also 
present us with exact expressions for the change in scale 
of cosmological solutions that has been predicted using 
Zalaletdinov's averaging scheme [5^H53] . More generally, 
we believe our study could allow some insight into the 
problem of how averaging should be performed in rela- 
tivistic cosmology. To date, most studies on this subject 
have used either highly symmetric exact solutions that 
contain a perfect fiuid, or have considered small fiuctu- 
ations around an FLRW "background" geometry. The 
former approach is strongly limited by the high degree 
of symmetry required in the solutions, while the latter is 
limited to geometries that are already necessarily close to 
FLRW. The models we have constructed here suffer from 
neither of these shortcomings, as they admit no Killing 
vectors, and do not require the assumption of an FLRW 



geometry, either as a boundary condition for the inhomo- 
geneities, or as a background geometry. As such it offers 
a new laboratory for testing ideas about inhomogeneity, 
averaging, and backreaction in cosmology. 

With the notion of effective mass employed here, we 
find the interaction energy grows rapidly as the number of 
masses in the lattice, and comes to dominate in the limit 
of very many masses. This means that if we compare 
lattices with different numbers of masses, but with the 
same total energy (including interaction energy), then 
the scale of the hypersurface of maximum expansion in- 
creases dramatically with increasing number of masses. 
As was discussed above, however, the interpretation of 
this effect requires some care, as it is based on a gener- 
alization of the definition of effective mass given by Brill 
and Lindquist in asymptotically flat settings. Neverthe- 
less, this effect is entirely neglected when treating the 
matter content of the universe as dust. 

Finally, in this paper, we have confined ourselves to 
the comparison of the discrete and continuous models 
on time symmetric hypersurfaces. Clearly the next 
step would be to make a detailed comparison of the 
full evolution of these models. We shall return to this 
question in future publications. We also note that, 
although we have only considered regular arrangements 
of masses in this paper, the formalism we have used 
allows for the possibility of considering much more 
complicated distributions. 
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Note added. - After submission of our manuscript the 
following related work appeared: [551 - 157] . The first of the 
these papers performs a numerical analysis of a spatially 
flat lattice of black holes [35] , while the second performs 
a perturbative analysis of a similar situation 36 . The 
third study finds exact initial data for a lattice of eight 
black holes in a space with spherical topology (using the 
same method presented in this paper) , and then proceeds 
to numerically evolve this data [37j . 
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